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An  automated  linear  measurement  sys- 
tem (LMS)  has  been  developed  to  de- 
termine the  nonlinearity  of  a tuned 
30  MHz  power  detector  over  a 6.021 
dB  range.  This  detector  uses  a single 
thermistor  bead  design  with  thermal 
isolation  to  obtain  nearly  linear  tracking 
over  a 4:1  change  in  input  power.  The 
nonlinear  correction  for  this  change, 
determined  by  the  LMS,  is  on  the  or- 
der of  1.00030  (-1-130  ^J,B)  for  the  de- 
tector presently  in  use.  Initial 


experiments  indicate  an  expanded  un- 
certainty of  ±0.138%  (±598  p-B), 
which  is  based  upon  Type  A and 
Type  B components. 
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1.  Introduction 

There  has  been  a recent  interest  in  and  demand 
for  a calibration  service  at  NIST  to  support  rf  at- 
tenuators and  voltage  doublers  that  operate  specif- 
ically at  30  MHz.  The  first  step  required  to  offer 
such  a service  is  to  develop  a reference  standard. 
For  the  best  possible  accuracy,  a tuned  single- 
element thermistor  mount  was  chosen.  A linear 
measurement  system  was  designed  and  constructed 
at  NIST  to  calibrate  the  nonlinearity  of  this  mount. 
This  paper  contains  a description  of  the  LMS,  an 
explanation  of  the  measurement  scheme,  calibra- 
tion results,  and  an  uncertainty  analysis. 

2.  System  Description 

A diagram  of  the  system  is  shown  in  Fig.  I.  The 
signal  generator  provides  a stable  30  MHz  rf  signal 
that  is  amplified  and  filtered  before  the  signal  is 
split  into  two  channels.  One  channel  consists  of  a 
variable  phase  shifter  and  a fixed  attenuator.  A 
coaxial  switch  either  terminates  the  signal  with  a 
50  n termination  or  feeds  the  signal  into  a power 
divider,  which  splits  the  signal  again.  Half  of  the 


signal  is  detected  by  a single-element  thermistor 
mount.  Pi,  and  the  other  half  is  fed  into  the  hybrid. 
The  second  channel  consists  of  a variable  attenua- 
tor. Like  the  first  channel,  a coaxial  switch  either 
terminates  the  signal  with  a 50  fi  load  or  feeds  the 
signal  into  a power  divider,  which  splits  the  signal. 
Half  of  the  signal  is  detected  by  a single-element 
thermistor  mount,  P2,  and  the  other  half  is  fed  into 
the  hybrid. 

The  hybrid  takes  the  sum  and  difference  of  the 
two  input  signals.  The  difference  is  fed  into  a diode 
detector  to  rectify  the  signal,  and  then  into  a null 
meter.  The  sum  is  fed  into  a coaxial  switch,  and 
when  the  switch  is  in  position  I,  the  signal  is  de- 
tected by  a third  thermistor  mount,  P3,  which  is  the 
thermistor  to  be  calibrated. 

Each  thermistor  mount  is  connected  to  a NIST 
Type  IV  bridge  and  a digital  voltmeter  to  measure 
rf  power. 

The  computer  controls  the  signal  generator,  the 
digital  voltmeters,  and  the  switch  controller,  and 
handles  the  data  acquisition  and  processing  through 
an  IEEE-488  bus. 
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Fig.  1.  Block  diagram  of  30  MHz  linear  measurement  system. 


3.  Design  of  30  MHz  Single-Element 
Thermistor  Mount 

The  dc  response  of  detector  P3  must  be  nearly 
linear  with  changes  in  input  rf  power.  A single  50  fl 
thermistor  bead  design  was  selected  as  a linear  de- 
tection scheme  [1].  This  detector  is  used  in  con- 
junction with  a NIST  Type  IV  self-balancing 
dc-substitution  rf  power  meter  modified  to  bias  a 
50  fi  detector  [2].  The  Type  IV  power  meter  is  de- 
signed to  change  the  thermistor  bead  bias  current 
so  that  the  thermistor  always  maintains  the  same 
resistance.  The  detectors  are  placed  in  the  LMS 
housing  where  the  temperature  is  held  to  ±0.2  °C. 

It  is  difficult  to  filter  the  rf  signal  from  the  power 
meter  leads  due  to  the  nature  of  the  single  bead. 
The  total  rf  signal  should  appear  across  the  bead  in 
the  ideal  case.  Other  considerations  in  the  detector 
design  include  thermal  stability  and  any  forms  of  rf 
leakage  into  or  out  of  the  detector. 

Several  ideas  are  incorporated  in  the  detector  to 
eliminate  these  problems.  Figure  2 shows  the  cir- 
cuit diagram  of  a single-element  thermistor  mount. 
LC  filter  sections  are  inserted  to  filter  the  rf  signal 


in  the  power  meter  leads.  No  ferrite  material  (such 
as  a ferrite  core  inductor)  is  used  near  the  thermis- 
tor, since  ferrite  components  experience  changes  in 
impedance  with  changes  in  the  rf  power.  These 
impedance  changes  lead  to  nonlinearities  in  the  de- 
tector response.  Therefore,  air-cored  inductors  are 
inserted  near  the  thermistor.  Ferrite-cored  conduc- 
tors are  allowable  in  sections  following  the  first  fil- 
ter section  because  the  rf  power  is  sufficiently 
reduced  and  renders  any  impedance  variation  neg- 
ligible. The  parallel  LC  filter  is  tuned  to  resonate 
at  30  MHz.  A special  thermistor-mounting  struc- 
ture has  been  developed,  and  consists  of  an  electri- 
cally insulated  doughnut  and  two  copper  blocks, 
one  on  either  side  of  the  doughnut.  The  thermistor 
is  placed  inside  the  doughnut  hole  and  is  encapsu- 
lated in  an  air  pocket  by  the  addition  of  the  copper 
blocks.  The  copper  blocks  also  provide  a large 
thermal  mass  so  that  the  entire  structure  cannot 
experience  rapid  changes  in  temperature.  The 
thermal  time  constant  is  much  longer  than  the  time 
required  to  perform  a single  measurement  cycle. 
Double-sided  copper-clad  fiberglass  boards  are 
used  in  the  exterior  detector  housing  and  in  the 


20 


Volume  99,  Number  1,  January-February  1994 

Journal  of  Research  of  the  National  Institute  of  Standards  and  Technology 


Power  Meter  Leads 


3 mW  will  be  incident  on  P3.  Finally,  switch  2 is 
moved  back  to  position  1 and  switch  1 is  moved  to 
position  2,  so  that  power  is  only  applied  to  the  sec- 
ond channel.  Readings  are  taken  on  P2  and  P3  and 
are  designated  Piou  and  P3oir,  respectively.  A nomi- 
nal 12  mW  will  be  incident  on  P2  and  approximately 
3 mW  will  be  incident  on  P3. 

The  calibration  constant  of  P3,  denoted  C hl,  is 
calculated  using 


Fig.  2.  Circuit  diagram  of  single-element  thermistor  mount. 
All  inductors  are  hand-wound  and  all  capacitors  are  1 jcF  unless 
otherwise  labeled. 


internal  compartment  walls  which  separate  various 
filter  sections.  This  type  of  construction  and  the 
use  of  capacitive  feedthroughs  in  the  compartment 
walls  reduce  the  amount  of  rf  leakage  though  the 
detector.  These  physical  and  electrical  construction 
considerations  provide  a stable  linear  detector  with 
low  rf  leakage. 

4.  Measurement  Methods 

4.1.  Calibration  of  P3 


The  derivation  of  this  formula  can  be  found  in 
Appendix  A.  The  calibration  constant  is  a measure 
of  the  detector’s  nonlinearity  over  the  6.021  dB 
power  change,  and  is  used  as  a multiplication  factor 
to  correct  the  ratio  measured  by  P3,  where 


2 Pi\U 
HL  -5 
i 3IOr 


(2) 


4.2.  Power  Measurements 

The  NIST  Type  IV  power  meter  must  be  con- 
nected to  an  external  dc  voltmeter.  The  substituted 
dc  power,  Pdc,  is  calculated  from  measured  voltages 
using 


Three  power  meters  are  used  in  the  LMS  — Pi,  P2, 
and  P3.  When  a given  power  meter  is  read,  the 
notation  used  is  Pxvzr,  where  X denotes  the  power 
meter  (1,  2,  or  3),  Y is  1 if  channel  1 is  switched  on 
and  0 if  channel  1 is  switched  off,  Z is  1 if  channel 
2 is  switched  on  and  0 if  channel  2 is  switched  off, 
and  “r”  denotes  that  this  power  is  a reading  and  not 
a true  value.  Powers  without  a subscript  “r”  are  true 
values. 

The  calibration  of  P3  is  achieved  using  a three- 
stage  method.  First,  with  switches  1,  2,  and  3 each 
in  position  1,  enough  power  is  applied  so  that  a 
nominal  12  mW  is  incident  on  the  three  thermistor 
mounts.  The  phase  shifter  and  variable  attenuator 
are  adjusted  to  balance  the  two  channels,  thus 
obtaining  a null  on  the  null  meter.  Readings  are 
taken  on  Pi,  P2,  and  P3  and  are  designated  Pnir, 
P2iir,  and  P3iir,  respectively.  Next,  switch  2 is  moved 
to  position  2,  so  that  power  is  only  applied  to  the 
first  channel.  Readings  are  taken  on  Pi  and  P3  and 
are  designated  Puor  and  P3ior,  respectively.  A nomi- 
nal 12  mW  will  be  incident  on  Pi  and  approximately 


“ On  ' on 

rdc-  B » (j) 

/\o 

where  Votf  is  the  output  voltage  with  no  rf  power  ap- 
plied, Von  is  the  output  voltage  with  rf  applied,  and 
Ro  (50  17)  is  the  resistance  of  the  thermistor  mount. 
Figure  3 shows  the  measurement  sequence  for  a 
power  calculation  [3].  An  initial  Koff  is  taken;  rf 
power  is  applied  and  Von  is  measured;  rf  power  is 
removed  and  a final  Votf  is  taken.  The  initial  and  fi- 
nal dc  measurements  are  used  with  the  Von  mea- 
surement to  calculate  the  power  and  correct  for  any 
mount  drift,  which  is  assumed  to  be  linear.  The  cal- 
culated value  of  Foff  in  Eq.  (3)  is  given  by 

Koff  =Koff,i  + ^ (Koff,f-l/of,i)  , (4) 

h — ‘i 

where  Koff.i  is  the  voltage  reading  taken  before  rf  is 
applied  at  time  ti,  Voff,f  is  the  voltage  taken  after  rf 
is  removed  at  time  (3,  and  t2  is  the  time  at  which  Kon 
is  taken. 
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TIME 


Fig,  3,  Sequence  for  measuring  power  meter  dc  voltages. 

5.  Results 

The  calibration  constant,  Chl,  which  is  a measure 
of  the  nonlinearity  of  mount  P3,  was  obtained  by  re- 
peated measurements.  The  average  value  of  the 
three  hundred  trials  taken  so  far  is  Chl  = 1-00030  or 
+ 130  |xB.  Long-term  data  are  being  accumulated 
to  validate  the  calibration  of  P3. 

Table  1 shows  a sample  calibration.  Powers  read 
by  the  three  thermistor  mounts  are  displayed  for 
each  of  the  three  stages.  The  calculated  calibration 
constant  is  shown,  along  with  the  actual  step  in 
power  for  each  leg  of  the  system. 


Table  1.  Sample  calibration  results  of  the  30  MHz  Linear 
Measurement  System’s  standard  mount 


Switches  1&2 

Switch  1 

Switch  2 

P(l) 

11.5194 

11.5089 

0.0000 

P(2) 

11.4170 

0.0001 

11.3991 

P(3) 

11.0907 

2.7644 

2.7756 

Leg  1: 

: 6.0335  dB 

Leg  2: 

: 6.0160  dB 

Chl  : 

1.00029 

6.  Uncertainty  Analysis 

6.1.  Evaluation  of  Type  A Standard  Uncertainty 

Evaluation  of  a Type  A standard  uncertainty  may 
be  based  on  any  valid  statistical  method  for  treating 
data.  Examples  are  calculating  the  standard  devia- 
tion of  the  mean  of  a series  of  independent  observa- 
tions, using  the  method  of  least  squares  to  fit  a 
curve  to  data  in  order  to  estimate  the  parameters  of 
the  curve  and  their  standard  deviations,  and  carry- 
ing out  an  analysis  of  variance  in  order  to  identify 
and  quantify  random  effects  in  certain  kinds  of 
measurements  [4]. 


The  calibration  of  the  standard  mount,  P3,  has 
been  repeated  three  hundred  times  to  determine 
the  repeatability  of  the  system.  Tests  were  per- 
formed at  various  times  of  the  day  over  several  days 
to  cover  as  many  random  factors  as  possible,  includ- 
ing variations  of  environmental  conditions  and  the 
operator’s  ability  to  renull  the  system.  The  sample 
standard  deviation  of  the  mean  is  ±0.000335%  or 
±2  |xB.  Long-term  data  are  being  accumulated  to 
validate  this  calibration,  and  a control  chart  is  being 
developed  to  track  any  possible  outliers  or  drift. 

6.2.  Evaluation  of  Type  B Standard  Uncertainty 

Evaluation  of  a Type  B standard  uncertainty  is 
based  upon  scientific  judgment  using  all  of  the 
available  relevant  information.  This  includes  previ- 
ous measurement  data,  manufacturers’  specifica- 
tions, data  provided  in  calibration  reports, 
knowledge  of  the  behavior  of  relevant  instruments 
and  materials,  and  uncertainties  assigned  to  refer- 
ence data  taken  from  handbooks  [4]. 

The  Type  B evaluation  of  standard  uncertainty 
accounts  for  the  following  factors: 

a.  Uncertainty  in  the  dc  voltmeter  measurements. 

b.  Uncertainty  in  the  Type  IV  power  meters. 

c.  Imperfect  isolation  between  the  two  legs. 

d.  Uncertainty  because /’301 5^ P310. 

e.  Effects  of  impedance  changes  in  P3. 

f.  RF  leakage. 

g.  Spurious  signals  and  harmonics. 

6.2.1  Voltmeter  Uncertainty  The  uncertainty 
in  the  individual  voltmeter  readings  may  be  deter- 
mined by  taking  the  total  differential  of  the  power 
expression,  Eq.  (3),  which  gives 

dP  = (Koff  dUoff-  Lon  dLon)  . (5) 

Ao 

The  total  differential  of  power,  Eq.  (5),  may  be  de- 
termined by  taking  the  differential  of  Koff,  Eq.  (4), 
which  gives 

dLoff  = (1  - r)dKoff,i  + rdKoff.f  , (6) 

where 

The  uncertainties,  dLoff.i,  dLoff,f,  and  dFon,  in  the 
measured  values  of  Loff.i,  Loff.f,  and  Lon,  are  based  on 
the  voltmeter  manufacturer’s  specifications. 

Figure  4 shows  the  uncertainty  in  the  power 
measurement  as  a function  of  power  level,  assum- 
ing the  powers  are  ratioed  as  they  are  in  the  Chl 
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calculation,  Eq.  (1).  The  power  measurements,  P30] 
and  P310,  which  are  approximately  3 mW,  result  in 
uncertainties  of  0.036%.  The  other  power  measure- 
ments, which  are  approximately  12  mW,  result  in 
uncertainties  of  0.008%.  The  uncertainty  of  Chl 
due  to  the  voltmeter  readings  may  be  found  by  in- 
serting the  individual  power  uncertainties  into  Eq. 
(1),  which  gives 


HL 


-l-  Av  — 


Pm  + ^310 

Pno  — ^110 

P 311  + 4I311 

/’lll-zllll 


1/2 


-H 


P 301  + Aj 

Pioi  ~ A 


1/2 


201 


Pill  + ^311 

P 211  ~ 4i211 


(8) 


This  results  in  an  uncertainty  of  Ay  = ±0.028%  or 
±122  jjlB. 


Fig.  4.  Power  measurement  uncertainty  due  to  DVM  when 
ratios  are  taken. 

6.2.2  Type  FV  Power-Meter  Uncertainty  The 

four  possible  sources  of  uncertainties  internal  to  the 
Type  IV  power  meter  are  the  reference  resistors, 
the  operational  amplifier  open-loop  gain,  input 
offset  voltage,  and  input  bias  current.  Larsen  has 
shown  that  the  uncertainties  due  to  the  Type  IV 
power  meters  are  negligible  compared  to  those  of 
the  voltmeters  [2]. 


6.2.3  Imperfect  Isolation  An  uncertainty  in 
the  calculation  of  Cm.  may  result  if  P|  and  P2  are  not 
zero  when  they  are  assumed  to  be.  The  thermistor 
mounts  are  square-law  detectors  and  are  not  suffi- 
ciently sensitive  to  determine  whether  Pi  and  P2  are 
low  enough  in  the  “off’  condition  to  avoid  signifi- 
cant errors.  This  uncertainty  is  derived  in  Appendbc 
B. 

The  corrected  formula  for  Cm.,  taking  imperfect 
isolation  into  account,  is 


+ 


Pmr 

PziiU 

PilU 

P21U 


1/2 


2 


1 

|l  + y 


where 


and 


^ QijAbi 

h 

Olio 


Ahi 

biiiiQn 


(9) 


(10) 

(11) 


Here,  yi  is  the  measure  of  Pi’s  contribution  to  Ciu. 
if  it  is  not  zero  when  it  is  assumed  to  be,  and  fii  is 
the  measure  of  P2’s  contribution  to  Chl  if  it  is  not 
zero  when  it  is  assumed  to  be.  The  terms,  bm  and 
6201  are  the  corresponding  voltages  to  the  powers 
Pno  andPaii,  respectively.  Assuming  Pi  10  and  P201  are 
12  mW,  6110  and  ^201  are  equal  to  0.7746  V.  The 
term,  Qn,  is  defined  in  Appendix  B and  its  value  is 
approximately  —1.  The  Ab’s  represent  the  b’s 
which  were  assumed  to  be  zero  in  Appendix  A.  If 
the  isolation  between  the  channels  is  65  dB,  as  is 
stated  in  the  power  divider  manufacturer’s  specifi- 
cations, then  Abi  = Ab 2 = 0.000436.  Using  nominal 
values  for  the  P’s  and  substituting  the  values  into 
Eq.  (9),  an  uncertainty  of  ±0.112%  or  ±488  |jlB  is 
obtained. 

6.2.4  Uncertainty  Because  P3oi^P3io  An  un- 
certainty in  the  calculation  of  Chl  may  result  if 
Pm^Pm.  The  derivation  of  this  uncertainty  can  be 
found  in  Appendix  C. 

The  measure  of  nonlinearity,  a,  of  thermistor 
mount  three  is 


1 ~ Chl 

Chl  Pu  — Pm 


(12) 
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where  Pu  is  the  reading  of  mount  three  at  the  low 
level  and  Pw,  is  the  power  reading  of  mount  three 
at  the  high  level.  Assuming  Chl=  1.0004  worst 
case,  Pu  = 2>  mW,  and  Ph,  = 12  mW,  a is  calculated 
to  be  3.334  x 10“^  The  ratio  KmIKm  is  given  by 


A!^301  _ 1 T aP 3()lr 

ATsio  1 + aP 3l0r 


(13) 


If  ^^30^  is  10%  greater  than  P3ior,  then 


Km  ^ l+a(3xl.l0) 
K-iw  1 + a(3) 


= 1.000010. 


(14) 


The  effect  of  this  being  nonunity  results  in  a cor- 
rected formula  for  Chl,  where 


= 1.000005  (15) 

using  nominal  values  for  the  P’s.  The  result  is  an  un- 
certainty of  ± 0.0005%  or  ± 3 [xB. 

6.2.5  Uncertainty  Due  to  Impedance  Changes  in 
Pi  An  uncertainty  in  the  calculation  of  Chl  may 
result  if  the  impedance  changes  in  P3  with  power 
level.  The  derivation  of  this  uncertainty  can  be 
found  in  Appendix  D. 

The  corrected  formula  for  Chl,  taking  impedance 
changes  in  P3  into  account,  is 


C 


2 _ 

HL  — 


(16) 


where  Pm  is  the  reflection  coefficient  of  mount 
three  at  the  high  level,  Pno  is  the  reflection  coeffi- 
cient of  mount  three  at  the  low  level,  and  Tlms  is  the 
reflection  coefficient  of  the  system  looking  into  port 
three. 

The  measured  value  of  Pno  was  actually  taken 
10  dB  below  Tsio  instead  of  the  6 dB  step,  so  the 
uncertainty  should  be  conservative.  Using  nominal 
values  for  the  P’s  and  the  measured  reflection  coef- 
ficients, the  calculated  uncertainty  is  ±0.014%  or 
± 62  |i,B. 


Impedance  changes  in  the  diode  detector  with  re- 
spect to  power  level  are  assumed  to  have  a negligi- 
ble effect  on  the  overall  uncertainty. 

6.2.6  Uncertainty  Due  to  RF  Leakage  It  is  dif- 
ficult to  assign  a quantitative  uncertainty  due  to  rf 
leakage.  In  order  to  reduce  leakage,  all  coaxial  ca- 
bles in  the  system  were  replaced  with  semirigid  lines 
and,  wherever  possible,  SMA  connectors  were  used. 
For  now,  the  uncertainty  due  to  rf  leakage  is  as- 
sumed to  be  ±10  p.B  or  ±0.0023%. 

6.2.7  Uncertainty  Due  to  Spurious  Signals  and 
Harmonics  After  the  signal  is  amplified  and  fil- 
tered, any  harmonics  are  at  least  —92  dBc,  while 
spurious  signals  are  no  greater  than  -84  dBc.  This 
will  result  in  a measurement  error  that  will  affect 
the  uncertainty  of  Chl. 

From  Eq.  (3),  the  calculated  dc  power  is 

Pdc  = 20  (F^ff  -K^„),  (17) 

where  Pdc  is  in  mW.  Assuming  Foff=l  V,  then 
Ton  = 0.63246  V for  a calculated  power  of  12  mW 
and  Von  = 0.92195  V for  a calculated  power  of  3 mW. 

A harmonic  at  - 92  dBc  results  in  an  uncertainty 
of  ±0.0025%  for  a voltage  measurement.  This  un- 
certainty translates  to  power  uncertainties  of 
(12  ±0.0004)  mW  and  (3  ±0.0008)  mW.  Using  the 
nominal  powers  to  calculate  Chl,  an  uncertainty  of 
±0.023%  or  ±100  |jlB  is  obtained. 

A spurious  signal  at  —112  dBc  results  in  an  un- 
certainty of  ± 0.00025%  for  a voltage  measurement. 
This  uncertainty  translates  to  power  uncertainties 
of  (12  ±0.00004)  mW  and  (3  ±0.000085)  mW.  Us- 
ing the  nominal  powers  to  calculate  Chl,  an  uncer- 
tainty of  ±0.0025%  or  ±11  p,B  is  obtained. 

6.2.8.  Overall  Type  B Uncertainty  For  each 
Type  B component,  an  estimated  range,  ±aj,  is 
given,  assuming  that  the  quantity  in  question  has  a 
100%  probability  of  lying  within  that  interval.  The 
quantity  is  treated  as  if  it  is  equally  probable  for  its 
value  to  lie  anywhere  within  the  interval.  Therefore, 
it  is  modeled  by  a rectangular  probability  distribu- 
tion. The  best  estimate  of  the  standard  deviation, 
Uj,  is 


(18) 


Table  2 shows  all  of  the  Type  B components  along 
with  their  corresponding  uncertainties  and  standard 
deviations.  The  overall  standard  deviation  of  the 
Type  B components,  calculated  using  the  root-sum- 
of-squares  method  (RSS),  is  ±0.069%  or  ±299  |xB. 
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Table  2.  Type  B eomponents  of  the  30  MHz  Linear  Measure- 
ment System  with  corresponding  uncertainty  ranges  and  stan- 
dard deviations 


Component 

Range 

(%) 

Standard  deviation 

(%) 

dc  voltage  measurements 

±0.028 

0.016 

Imperfect  isolation 

±0.112 

0.065 

Aoi 

±0.0005 

0.0003 

Impedance  changes  in  P3 

±0.014 

0.008 

rf  leakage 

±0.0023 

0.0013 

Spurious  signals 

±0.0025 

0.0014 

Harmonics 

Combined  Type  B Standard 

±0.023 

0.013 

Uncertainty  (RSS) 

± 0.069 
(±299  pB) 

6.3  Combined  Standard  Uncertainty 

The  combined  standard  uncertainty,  Uc,  is  taken 
to  represent  the  estimated  standard  deviation  of  the 
result.  It  is  obtained  by  combining  the  individual 
standard  deviations,  u\,  whether  arising  from  a Type 
A or  a Type  B evaluation  [4].  The  technique  used  to 
combine  the  standard  deviations  is  the  RSS 
method. 

The  total  uncertainty  reported  is  the  expanded 
uncertainty,  U , which  is  obtained  by  multiplying  Uc 
by  a coverage  factor.  A:.  To  be  consistent  with  cur- 
rent international  practice,  the  value  of  k used  at 
NIST  for  calculating  U isk=2.  The  total  expanded 
uncertainty,  U,  of  the  mount’s  nonlinearity  is  calcu- 
lated to  be  ±0.138%  or  ±598  [xB. 

7.  Conclusion 

The  first  step  toward  offering  a calibration  ser- 
vice at  NIST  to  support  rf  attenuators  and  voltage 
doublers  that  operate  at  30  MHz  has  been  com- 
pleted. The  reference  standard,  a tuned  single- 
element thermistor  mount,  has  been  calibrated  us- 
ing an  LMS,  designed  and  constructed  at  NIST. 
The  next  step  is  to  modify  the  LMS  so  that  a device 
under  test  may  be  inserted  into  the  system  and  cal- 
ibrated against  the  reference  standard. 

8.  Appendix  A.  Scattering  Coefficient 
Analysis  of  the  LMS 

The  30  MHz  Linear  Measurement  System  can  be 
considered  a five  port  network,  as  shown  in  Fig.  5. 
The  analysis  assumes  that  the  LMS  is  an  ideal  sys- 
tem—that  is,  the  system  is  linear,  and  when  bi  and 
bi  are  assumed  to  be  0,  they  truly  are. 

The  pertinent  scattering  coefficients  are 


f)l=5iini+5'i2<32  + 5'i3^r3  + 5i4<34  + 5i5fl5 
^2  = 521^  1 + 522^^  2 + SiyCl  3 + 524^  4 + S25CI  5 
^3  = ‘^31^  1 + 532^1  2 + 5^33^  3 + ‘^34^  4 + ^35^^  5 

64  = 5'4lfl  1 + 2 ± 3 + SmO  4 ± ‘545^  5 

b5  = S5\Cl  1 +552fl2  + 5'53n3  + 5'54^4  + ‘S'5542  5 • (19) 

If  = Fifei  and  ai  = r2b2,  then  Eq.  (19)  may  be  writ- 
ten as 

bi(l  — T15ii)  = 512/2^2  + 5i3fl3  + ‘S'i4a4  + 

^2(1  “ T2522)  = ^21  Fibi  +523^3  + *52404  + 52505 
f>3  = 531  Fib  I +532^2^2  + 53303  + 53404  + 53505 
^34  = 541  Fibi  + 542T2  62  + 54303  + 54404  + 54505 

65  = 551  F161  + 552F262  + 553O3  + 554O4  + 555O5 

(20) 

This  gives  five  equations  and  eight  unknown  vari- 
ables—03,  04,  Os,  b\,  62,  63,  64,  and  65.  Four  variables 
can  be  eliminated  from  Eq.  (20),  say  04,  64,  05,  and 
65.  Eq.  (20)  becomes 

b\=P'\^i  + Q[2b2-^Q[ib2  , (21) 

where  the  primed  quantities  are  functions  of  the 
scattering  coefficients  and  Fi  and  r2. 

It  is  assumed  that  Fs  is  constant  with  respect  to 
varying  power  levels.  If  now  <33  = F363, 

6i  = (Fi'3  F3  + (2 13)63  + C 1262  = 

C1363  + C1262  (22) 

or 

b,  =*1^^ , (23) 

L:I3 

where  two  complex  constants,  Q12  and  Qn,  describe 
the  relationship  among  61,  62,  and  63. 


Fig.  5.  LMS  port  nomenclature. 
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Equation  (23)  may  be  applied  under  three  condi- 
tions: subscript  “10”  when  ^>2  = 0;  “01”  when  bi  = 0; 
and  “11”  when  neither  bi  nor  bz  equal  0.  Thus 

I biii  — Qi2b2n 

biu= , (24) 

where  £>311  is  the  value  of  £>3  when  bi  and  £>2  are  both 
energized  from  ports  four  and  five.  When  £>2  = 0, 
which  occurs  by  energizing  port  four  and  terminat- 
ing port  five,  Eq.  (23)  becomes 

bm=  ^ (25) 

Wfl3 

or 

^ • (26) 

D310 


Similarly,  when  £>i  = 0,  which  occurs  by  energiz- 
ing port  five  and  terminating  port  four,  Eq.  (23) 
becomes 

£>301  = “7^^  ^201  • (27) 

W'13 


Substituting  Eq.  (26)  into  Eq.  (27)  gives 


or 


6301  — 


-Qnb 


— Qi2  — 


bioi  biw 
£>201  £>310 


(28) 

(29) 


Substituting  Eqs.  (26)  and  (29)  into  Eq.  (24)  gives 


£>111  £>310  ^ £>211  £>301 
£>311  £>110  £>311  £>201 


(30) 


The  b’s  are  complex  numbers.  If  the  phase  is  mea- 
sured relative  to  £>iii/£>3ii,  then 


Since  the  detectors  measure  only  magnitudes,  the 
magnitude  squared  of  Eq.  (32)  is 


2 

sin^  6 , (33) 


which  can  also  be  written  as 


£>301 

£>201 

£>311 


£>310 

2 

£>301 

£>110 

+ 

£>201 

£>311 

£>311 

£>111 

£>211 

£>310 

£>301 

£>110 

£>201 

£>311 

£>311 

£>111 

£>211 

(34) 


Since  the  powers  are  proportional  to  the  squares  of 
the  voltages,  Eq.  (34)  becomes 


£>310 

£>110 


eJ»im  + 


Completing  the  square  and  transposing  terms  gives 


£>211 

£>311 


gj  ®2311 


£>301 

£>201 


q]  ^201 


(31) 


U 6 = 02311  -I-  03201  - 03110,  Eq.  (31)  becomes 


l=e^  03110 


’£>111 

£>310 

.£>311 

£>110 

£>211 

£>301 

£>311 

£>201 

(32) 
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Since  the  powers  are  all  positive  and  (1  — cos0)  is 
greater  than  or  equal  to  0,  the  quantity  on  the  left 
side  of  Eq.  (36)  is  a minimum  when  a null  is  ob- 
tained (0  = 0).  Since  cos0  approaches  1 as  0 ap- 
proaches 0,  Eq.  (36)  becomes 


where  Chl  is  the  calibration  constant  that  one  must 
multiply  the  PiwJPiw,  reading  by  to  get  the  true 
ratio 


bn 

Pin 


(42) 


1 = 


The  calibration  constant  should  be  used  only  when 
P3  is  operating  at  the  two  levels  for  which  it  was 
calibrated. 


The  above  P’s  are  the  true  powers.  The  detec- 
tors, however,  give  readings  proportional  to  the 
true  values.  Furthermore,  P3  could  have  a different 
proportionality  constant  for  the  two  levels  if  the 
detector  is  nonlinear.  The  ratio  of  these  two  pro- 
portionality constants  is  the  desired  result. 

If  the  subscript  “r”  represents  proportional  read- 
ings, and  the  nonsubscripted  P’s  represent  the  true 
values,  then 

Pi  = A:iPir 
P2  = A:2P2r 

P 311  =/C3ii  P3iir 

P310  — /C31oP310r 

P 301  = ^301  P 301r  (38) 


9.  Appendix  B.  Imperfect  Isolation 

It  is  desirable  to  have  an  estimate  of  the  error 
when  Pi  and  P2  are  not  0 when  they  are  assumed  to 
be.  The  detectors  are  square-law  detectors  and  are 
not  sufficiently  sensitive  to  determine  whether  Pi 
and  P2  are  low  enough  in  the  “off’  condition  to 
avoid  significant  errors. 

Let  Ab's  represent  the  b's  which  were  assumed 
to  be  0 in  Appendix  A.  Equation  (24)  is  still  valid, 
but  Eq.  (25)  should  now  be 

L bim  — QnAbi  bm  ( ^ Qn  Abi  \ 

— W, — “ & 1 

and  Eq.  (27)  should  now  be 


The  powers,  P310  and  P301,  will  ordinarily  be  very 
close,  and  essentially 


Ab\  — ^12  bun  _ _ 12126201  / , Ab\  \ 
Q\i  Q\i  \ 6201  !2i2  / 


K^m  — Km  (39) 

since  the  detectors  are  assumed  to  be  very  nearly 
linear.  Pi  and  P2  remain  essentially  constant.  Sub- 
stituting Eqs.  (38)  and  (39)  into  Eq.  (37)  gives 


1 = 


PsiOr 

PllOr 

P31lr 

Plllr 


-I- 


( 


P 301r 

P201r 

P 311r 

P211r 


(40) 


or 


/-2  -^^311 

Chi  = 


Km 


{ 


PilOr 
P llOr 

PiW! 

Pnu 


PiOU 

PltiU 

PjWt 

Pl\u 


(41) 


Let 


and 


r,  Qi2  Ab2 

Pi —T 

Olio 


= 

620.^12  • 


Combining  Eqs.  (43)-(46)  gives 

Oij  = N(l  + ft) 

O310 


and 


..  6301  6110  1 + 13\ 
6201  6310  1 -f-  yi 


Eq.  (24)  now  becomes 


1 = 


6111  6211  6301  6110  \ + (ij 

bill  6311  6201  6310  1 -H  yi 

^(i+w 

O310 


(44) 

(45) 

(46) 


(47) 


(48) 


(49) 
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or 


J _ ^111  ^310  1 ^ ^>211  ^301  1 


^311  ^110  1 + i3i  b^w  bioi  1 + yi 


(50) 


The  b's  and  /3i  and  y\  are  complex  numbers.  If  the 
phase  is  measured  relative  to  bm/bsu,  then 


1 = 


bill 


-'311 


^310 

bllO 


a;(93110-  Sp\) 


ll+^. 


+ 


^211 

b3ii 


82311 


Completing  the  square  of  Eq.  (56)  and  transposing 
terms  gives 


P 310 
Elio 

Em 


1/2 


1/2 


|l  + /3i 


+ 


1 


Pm 

Pun 

E311  I |l  + Ti| 
E211 


1+E«,  [l-cos(0  + 0p,-0y,)] 


(58) 


6301 

g><83201-8,l-J 

(51) 

bm 

|l  + r,|  ' 

where 

and 

— 

Arg  (l  + /3i) 

(52) 

0yi=^^  (1  + yi)  • 

(53) 

Equation  (51)  may  also  be  written  as 

6310 


2 — g/(%iio-e^i) 


'110 


•-  Oiii 


+ 


»/(  8+  8^1  - 8yj) 


^311 

6211 


|i  + yi|  -* 


where 


0 — 02311  + 03201  — 0311 


(54) 


(55) 


As  in  Appendix  A,  taking  the  square  of  the  magni- 
tude and  changing  to  power  gives 


1=- 


E310 

Elio 


^il  + /3i 


-f 


E301 

E201 


-I- 


|i  + yip 


Em  E211 

Esq  COS(0  + 0^1-  0y,), 


(56) 


where 


As  in  Appendix  A,  using  the  argument  on  the  func- 
tion being  minimum  when  the  angle  is  0,  and  sub- 
stituting in  the  proportional  readings 


(^2  _ 

HL  — 


1/2 


r llOr 

E311r 

E iiir 


|l  + i3i 


-I- 


/ E301r 

y/2  -| 

1 E201r 

1 1 

\ Ejllr 

/ |i  + yi| 

E2llr 

/ J 

(59) 


10.  Appendix  C.  Error  Due  to  P^oi^Pno 

Since  the  detectors  are  very  close  to  being  pre- 
cisely linear,  it  is  reasonable  to  model  the  nonlin- 
earity by 

E=AEr(l+«Er), 

where  a is  a measure  of  the  nonlinearity.  The  non- 
subscripted  E is  the  true  value  and  Er  is  the  propor- 
tional reading.  Substituting  Eq.  (60)  into  Eq.  (42) 

Eh  Enr  KPut  + (xPHr)  //:i\ 

K-‘^'"-K:-KPuO  + aPu) 


or 


,2  1 + «EHr 


Chl  — 


1 -I-  aPu 


Solving  for  a 


1-C^ 


HL 


C hl  Pu~P\\i 


(62) 


(63) 


Esq=  2 


E310  E301 
Elio  E201 


^ll-fEiPIl  + yiP 

mil  mil 


(57) 


The  ratio  KmlK^w  may  be  found  using  Eq.  (60) 


E301  _ K301P 30ir  _ KP-jau  (1  + aP 3oir) 

E310  KiwPiOU  /fEsiOr  (1  + ttEsiOr) 


(64) 
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K301  1 + OcP 301r 

K310  1 + OiP 3i0r 


(65) 


The  effect  of  this  non-unity  is  obtained  by  substi- 
tuting Eq.  (38)  and  (65)  into  Eq.  (37),  which  gives 


bi  — Q12  62  + ( — /lms  Q\i  1 3 + Qh)  bi  (70) 
or 

63(l-rLMsr3)=p^-^  . (71) 

!J13  ^13 

As  before,  making  b\  and  62  successively  0,  gives 


(2.3  = 


611 


6310  1 — /lMS/310 


and 


(72) 


a:3oi  p 


301  r 301r 
^310  P 201r 

PlUr 

T’zilr 


<2 12 
(2.3 


6301  / -| 

— T ft  —2  LMS 

6201 


Pm). 


(73) 


(66) 


Assuming  PxQ  = rm,  since  6310=6301,  gives 


or 


C 


2 _ 

HL  — 


/C311 

^310 


-^  = ^(l-TLMsr3io).  (74) 

fjn  0201 

Substituting  Eqs.  (72)  and  (74)  into  Eq.  (71),  and 
rearranging  terms,  gives 


^301  Pimi 
K-3W  P20U 

7^31 Ir 

P2\U 


(67) 


1 — 71ms  All  _6iii  6310  I 6211  6301 
1 — 7lms73io  63116110  bjii  b2m 


(75) 


Since  the  right  side  of  Eq.  (75)  is  identical  to  the 
right  side  of  Eq.  (30),  the  same  manipulations  may 
be  performed  to  yield 


11.  Appendix  D.  Effects  of  Impedance 
Changes  in  P3 

Starting  with  Eq.  (21), 


61 —7^1303  + (2 12  62  + <2 1363.  (68) 

If  ports  four  and  five  are  terminated  in  impedances 
that  null  the  hybrid,  61  and  62  are  almost  0 even 
when  fl3  is  energized  since  the  two  legs  of  each 
power  divider  are  isolated.  Thus 


1 — Elms  E311 

J 

6310 

bill) 

6301 

6201 

1—7  LMS  7 310 

\ 

6311 

6311 

\ 

6111 

6211 

6310 

6301 

1 ® 

6201 

6311 

6311 

6111 

6211 

(cos  0 - 1)  . 


(76) 


ri,Ms=^=-^,  (69) 

^3  (2.3 

where  Elms  is  the  reflection  coefficient  looking  into 
port  three  when  ports  four  and  five  are  terminated 
with  matching  impedances.  The  second  part  of  Eq. 
(69)  is  obtained  from  setting  61  and  62  to  0 in  Eq. 
(68). 

Substituting  Eq.  (69)  into  Eq.  (68),  and  using  03 
= Pi  63  yields 


As  before,  the  powers  are  proportional  to  the  6^’s 
except  for  P3,  because  of  the  change  in  p.  Thus 

E3ii  = L%^(l-|r3.iP)  , (77) 


E310  = I-^(1-|E31oP),  (78) 

^0 
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and 

P30.=-L^(i-|r30,P) 

= L^(i-|rMop),  (79) 


type;  J.  Juroshek  and  P.  Leyva  for  reviewing  this 
paper;  P.  Voris  and  K.  Talley  who  helped  assemble 
the  LMS;  G.  Rebuldela  for  his  advice  throughout 
the  project;  and  N.  Larsen  for  supplying  the  data 
found  in  Fig.  4. 
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P310 
Pi  10 
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) 


1/2 


-1- 
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